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-  ACT 

Incoherent  infrared  scattering  of  silicon  is  reviewed  and  scat- 
tering from  a  crystal  grown  along  a  <^110^>- axis  has  been  examined.   It 
is  concluded  that  the  incoherent  infrared  scattering  arises  from  regions 
in  the  cyrstal  bounded  by  dislocation  loops.   The  observed  radii  of 
cyration  of  the  scattering  regions,  —  5  ix   ,  agree  with  the  size  of 
dislocation  loops  which  are  observed  using  Dash's  decorating  technique. 
Analysis  of  the  small  angle  scattering  precludes  the  presence  of  oxide 
particles  smaller  than  about  IjU  ,    and  microscopic  examination  discloses 
no  particles  larger  than  this  indicating  that  there  has  been  no  coalesence 
of  oxides. 
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CHAPTER  I 

INTRODUCTION 

Infrared  radiation  has  been  employed  by  many  investigators  in  at- 
tempts to  study  imperfections  in  the  structure  of  silicon  crystals. 
The  majority  of  efforts  have  been  directed  toward  the  deliniation  of 
dislocation  and  impurity  defects.   Chemical  etching  has  long  been  used 
as  a  means  of  visual  detection  of  the  emergence  of  dislocations  at 
crystal  surfaces.   A  method  of  creating  images  of  dislocations  within 

the  bulk  of  the  crystal  by  a  decorating  technique  has  been  developed  by 

1  2 
Dash  '   at  the  General  Electric  Laboratories.   He  was  able  to  trace 

visually  the  dislocation  lines  through  the  crystal  and  detect  their  e- 

mergence  at  etch  pits  at  the  surface.   X-ray  extinction  microscopy  has 

3  4 
been  applied  in  determining  dislocation  densities.  '    A  unique  method 

using  an  integrating  sphere  to  detect  light  scattering  by  both  disloca- 
tions and  oxide  impurities  was  used  at  the  General  Telephone  and  Elec- 

4 
tronics  Laboratories. 

This  investigation  was  directed  toward  the  application  of  incoherent 
scattering  techniques  as  a  means  of  particle  size  determination  in  sili- 
con crystals.   Small  angle  X-ray  diffraction  is  a  well  established 
method  used  to  measure  particle  sizes  of  the  order  of  magnitude  of  X-ray 
wavelength.   However,  this  method  is  not  applicable  to  structure  deter- 
minations of  sizes  large  compared  to  the  wavelength  of  the  radiation  used. 
Theoretically,  with  longer  wavelengths,  one  could  determine  the  sizes  of 
larger  particles.   The  problem  encountered  is  that  softer  X-rays  and 
visible  light  are  attenuated  to  quickly  by  matter.   The  wavelengths  in  the 


near  infrared  spectrum  are  of  the  magnitude  of  dislocation  and  impurity 
networks  observed  by  the  methods  referred  to  above.   Therefore,  in- 
coherent scattering  of  infrared  appears  to  offer  possibilities  as  an- 
other means  of  determining  the  dislocation  densities  within  silicon 
crystals.   It  is  also  proposed  that  a  diffraction  system  can  be  built  easily 
and  cheaply  thus  avoiding  the  highly  designed  and  costly  equipment  associ- 
ated with  X-ray  diffraction  and  small  angle  cameras. 


en  .   ii 

REVIEW  OF  DISLOCATION  AND  OXIDE  IMPURITY  DETECTION  IN 

SILICON  CRYSTALS 

1.  Chemical  Etching  -   A  method  of  identifying  dislocations  in  LiF 

5 
crystals  by  observations  of  etch-pits  reported  by  Gilman  has  been 

applied  to  semi-conductor  materials  by  other  investigators  (Vogel  et  al  , 

7  8 
Dash  '  ).   Dash's  work  was  devoted  to  the  study  of  dislocations  in  sili- 
con crystals.   He  reported  etch-pits  several  hundred  microns  in  diameter 
extending  a  comparable  distance  into  the  crystal  when  etched  overnight 
with  an  echant  consisting  of  the  following  parts  by  volume:   1  50%  HF, 
3  70%  HNO  ,  10  glacial  acetic  acid.   One  can  obtain  the  dislocation 
density  of  a  crystal  by  counting  the  etch-pits  which  mark  the  emergence 
of  dislocations  at  the  crystal  surface.   Dash  observed  deep  pits  on  any 
surface  independent  of  its  crystal lographic  direction.   The  disloca- 
tion density  within  any  one  crystal  may  vary  strongly  with  orientation, 

and  etch-pit  counts  on  only  one  surface  can  lead  to  erroneous  conclusions 

4 
about  crystal  perfection.   Schwuttke  reported  errors  up  to  100  per  cent 

in  crystals  of  low  dislocation  density.   Etch-pit  counting  is  universally 

applied  by  investigators  as  a  means  of  obtaining  dislocation  densities  in 

silicon  crystals. 

7  8 

2.  Decoration  -  Dash  *   describes  a  technique  for  observing  visually 

the  dislocation  lineage  within  the  bulk  of  silicon  crystals.   Copper 
atoms  diffused  into  the  crystal  were  precipitated  in  aggregates  along 
dislocations  after  appropriate  heat  treatment.   Silicon  but  not  copper, 
is  transparent  to  neai  infrared  radiation.  Using  an  infrared  image  tube 


fitted  with  a  microscope  he  was  able  to  obtain  images  of  dislocation 
lineages  in  the  crystals.   A  one  to  one  correlation  between  the  images 
and  dislocations  was  obtained  by  tracing  the  lines  to  etch-pits  at 
crystal  surfaces.   Precipitates  are  frequently  not  found  on  pure  screw 
dislocations,  but  when  formed  they  appear  as  tightly  wound  helicies. 
In   <3.1(£>  grown  crystals  he  observed  dislocation  loops  ranging  up  to 
50  microns  in  diameter.   Copper  atoms  appear  as  a  fine,  speckeled  sub- 
structure in  dislocation- free  crystals.  This  is  apparently  a  result  of 
agglomeration  of  the  precipatate  throughout  the  crystal  body. 

3.   Light  Scattering 

4 
a.   Integrating  Sphere  Technique  -  Schwuttke   reports  a  method  of 

detecting  incoherent  scattering  of  infrared  by  "bad  regions"  within 
silicon  crystals.   He  calculated  the  number  of  scattering  centers  (re- 
gions outlined  by  dislocations)  necessary  to  produce  detectible  scatter- 
ing.  Three  assumptions  were  made:  1.   Imperfection  scattering,  to  be 
observable,  must  be  ten  t.imes  greater  than  thermal  scattering,   2.   The 
scattering  centers  are  spherical  in  shape,  3.  The  radius  of  the  scatter- 
ing centers  is  of  the  order  of  0.1  micron.   Since  the  particle  size  is 
small  comparted  to  the  wavelength  of  the  radiation  (  }\   =  1. 1  JUL   ),  ex- 
pressions for  Rayleigh  scattering  are  developed  relating  intensity  to 

the  size  and  number  of  bad  regions  per  unit  volume.   Theoretical  consider- 

8       2 
ations  suggest  that  a  dislocation  density  of  at  least  10  per  cm  is 

necessary  to  produce  observable  scattering.   This  was  confirmed  in 
Schwuttke* s  experiment,  but  the  results  do  not  constitute  a  confirma- 
tion of  the  assumed  particle  size.   The  object  of  the  study  was  to 
determine  the  dislocation  density  necessary  to  produce  scattering. 


As  previously  noted,  Dash  reported  dislocation  loops  of  from  2-50  microns 
in  diameter.   Also  the  loops  tend  to  lie  in  favorably  orientated  planes 
thus  suggesting  a  platelet  configuration.   Regions  of  0.1  micron   (the 
assumed  particle  size  in  this  study)  described  by  dislocations  cannot  be 
resolved  by  Dash's  method,  but  since  the  scattering  depends  on  a  radius 
of  gyration  there  exists  an  infinite  number  of  particle  size  and  shape 
combinations  which  could  produce  the  same  scattering  effect.   Lacking 
supplementary  information,  one  cannot  draw  definite  conclusions  about 
either  size  or  shape. 

Schwuttke's  ingenious  experimental  technique  of  observing  very  small 
intensities  of  scattered  radiation  warrants  a  brief  description.   The 
optical  arrangement  of  the  instrument  appears  in  figure  2-1, 

A  gold-plated  integrating  sphere  contains  a  tungsten  filiment  light 
source  S.   The  black  body  B  and  openings  0  and  09  define  a  conical  umbra 
into  which  no  light  originating  inside  the  sphere  can  penetrate.   How- 
ever, scattering  by  a  crystal  placed  outside  of  0~  will  permit  light  to 
enter  the  "dark"  area.   With  a  suitable  optical  system  and  infrared  image 
tube,  scattered  light  can  be  observed. 

b.   X-ray  Extinction  -  X-ray  diffraction  techniques  have  been  used 
to  substantiate  dislocation  counts  in  crystals  obtained  by  the  previous- 
ly described  techniques.   The  double-crystal  spectrometer  rocking  curve 

9 
method  is  commonly  applied  to  determine  crystal  perfection.   Lang   re- 
ports this  method  to  be  unreliable  for  dislocation  densities  below  about 
10  per  cm  ,   X-ray  topographs  can  present  an  unambiguous  display  of 
dislocation  lineage  within  a  crystal.   He  was  able  to  observe  directly 
individual  dislocations  by  copper  atom  decoration.   Lang's  results  sub- 
stantiate the  preference  of  the  \HQx  direction  for  straight  line 
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Fig,  2-1   Schwuttke's  Scattering  Integrating  Sphere 


dislocations  reported  by  Dash.   Arrays  of  hexagonal  loops  are  formed 
by  the  closing  of  \11(£>  edge  dislocation  lines  by  pure  screw  component s . 
This  suggests  a  thin  two  dimensional  platelet  type  of  configuration  for 
a  particle  shape,         • 

c.   Additional  Results  -  One  ircoinplete  report  of  dislocation 
observation  was  found  in  the  literature.   The  results,  but  apparently 
not  the  experimental  procedures,  were  published.   Patel,  Wagner  and  Moss 
have  observed  the  aggregation  and  collapsing  of  supersaturated  vacancies 
to  form  condensed  prismatic  dislocation  loops  upon  cooling  silicon  crystals 
from  the  melt.   They  report  that  indirect  evidence  of  this  was  obtained 
by  Chandhuri  and  Patel  from  plastic  deformation  experiments.   The  dis- 
location loops  observed  after  deformation  were  estimated  to  be  about  1 
micron  in  diameter. 

4.   Oxygen  Impurities  -   Oxygen  impurities  have  been  observed  by  both 
decoration  and  light  scattering  techniques.   Several  investigators  have 
reported  the  effect  of  oxygen  in  contributing  to  erroneous  dislocation 
etch-pit  counts  but  have  not  used  echants  as  a  means  of  oxide  detection. 

Debye  and  Lederhandler   in  an  incomplete  report  state  they  have 
obtained  SiO~  particle  sizes  by  measuring  angular  dependence  of  scatter- 
ed infrared  radiation.   They  found  particle  sizes  of  the  order  of  the 

wavelength  of  the  light.   The  sizes  are  strongly  dependent  upon  the  heat 

12 
treatment  given  the  specimens.   Kaiser   reports  the  effect  of  long  (50 

hours)  heat  treatment  at  1000°C  on  the  optical  properties  of  silicon. 

i : 
Prior  to  heating,  silicon  exhibits  a  sharp  infrared  absorption  band  at 

a  wavelength  of  9  microns.   The  absorption  is  proportional  to  oxygen  con- 
centration.  Following  heat  treatment,  the  band  decreases  but  without  an 


accompanying  lowering  in  :ation.   He  attributes  this  to 

the  clustering  of  silicon  oxide  into  particles  of  lower  dielectric 
constant  than  the  silicon  matrix. 

Kaiser  also  makes  an  attempt  at  SiCL  particle  size  determination 
by  observation  of  Rayleigh  scattering.   He  assumes  a  spherical  particle 
of  size  smaller  than  the  wavelength  of  radiation  employed.   Since  Si  and 
SiO~  do  not  absorb  at  wavelengths  of  l-3yU  ,  changes  in  intensity  will 
be  due  to  Rayleigh  scattering  and  not  absorption.   For  spherical  parti- 
cles, the  absorption  coefficient  (  oC  )  due  to  Rayleigh  scattering  is 

oC  -  KMp  vXi4"  2-1 

where  K  is  a  constant,  Np  the  number  of  scattering  centers,  V  the  total 
volume  of  scattering  centers  and  A  the  wavelength.   Following  the  in- 
verse 4th  power  lav/,  a  log-log  plot  of  oC  vs.  A  should  yield  a  straight 
line  of  slope  equal  to  -4.   Experimental  results  produced  a  slope  of  -3.4 
which  adequately  substantiates  the  existance  of  Rayleigh  scattering. 
From  2-1  and  the  auxiliary  equation 

Nc 


No  /  M  \ 


3 
where  N   is  the  total  number  of  oxygen  atoms  per  cm  ,  M  the  molecular 

wright  of  SiO?,  ^  the  density,  and  A  avogradro's  number,  the  number  of 

scattering  centers,  Np  can  be  calculated.   The  number  of  oxygen  atoms 

per  cluster,  No/Np,  is  obtained  and  from  the  assumption  of  spherical 

particles  he  calculates  the  diameter  of  scattering  to  be  0.1  M.   This 

value  is  consistant  with  the  magnitudes  observed  by  Debye  and  Lederhand- 

ler  from  angular  dependence  scattering  experiments. 


Schwul  >st  conclusive  and  significant 

information  about  the  distribution  of  oxygen  in  silicone   He  calculated 

10  3 

that  a  minimum  of  10   oxygen  atoms  per  cm  was  necessary  to  produce 

scattering.   This  was  confirmed  by  observations  with  his  integrating 
sphere  technique.   An  investigation  was  .then  conducted  to  determine 
the  distribution  of  the  oxide  in  the  crystal.   By  application  of  etch-pit 
counts,  decoration  and  light  scattering  techniques  to  the  same  specimens 
he  was  able  to  make  observations  of  the  effect  of  dislocations  upon  the 
oxide  particles.   It  was  found  that  in  decoration  by  copper  atoms,  the 
precipitation  is  influenced  by  the  oxygen  content  of  the  crystal.   This 
suggests  that  oxygen  is  present  in  the  dislocation  networks  and  acts  as 
nucleation  centers  for  the  copper  precipitates.   The  copper  prefers  re- 
gions of  high  oxygen  content  thus  delineating  high  density  oxide  clusters, 
Dislocations  offer  favorable  sites  for  Si0?  clusters.   In  relatively  dis- 
location free  crystals,  clustering  by  heat  treatment  tends  to  take  place 
in  the  bulk  of  the  crystal, 

5.   Summary  -  From  the  foregoing  review  the  following  conclusions  are 
presented: 

1.  Dislocation  networds  will  produce  incoherent  infrared  scatter- 
ing and  allow  one  to  make  qualitative  conclusions  about  their  distribu- 
tion and  size. 

2.  The  assumption  of  a  spherical  shaped  particle  size  for  dis- 
location regions  is  not  necessarily  valid.   Decorating  techniques  suggest 
a  fishnet  structure  with  a  particle  shape  approaching  thin  platelets. 

3.  Particles  formed  by  dislocation  networks  have  some  degree  of 
preferred  orientation. 


4.  Oxygei  •        concentrate  on  dislocation  sices.   H< 
treatment  at  1000°C  produces  higher  concentrations  of  SiO- 

sites  through  clustering.   This  constitutes  an  increase  in  particle  size, 

5.  In  crystals  of  low  dislocation  density,  SiO-  molecules  will 
cluster  in  the  body  of  the  crystal.   Assumptions  of  a  sjherical  cluster 
shape  are  not  substantiated. 

6.  Dislocations  and  oxide  impurities  of  appropriate  size  are 
present  to  cause  incoherent  scattering, 

7.  Direct  evidence  exists  of  dislocation  loops  as  large  as  50yU 
in  diameter.  Strong  indications  exist  for  sizes  as  small  as  1AA  dia- 
meter. 


10 


CHAPTER  III 

THEORY  OF  SMALL-ANGLE  SCATTERING  OF  X-RAYS 

A,   General  -  The  purpose  of  this  study  is  to  investigate  the  feasi- 
bility of  applying  small-angle  scattering  techniques  to  the  determina- 
tion of  dislocation  networks  or  impurity  concentrations  in  silicon 
crystals,   These  heterogeneities  constitute  a  particle  distribution 
within  the  material.   Therefore,  it  is  appropriate  that  a  review  of  the 
theory  of  small-angle  scattering  by  particles  be  presented  at  this  point, 

The  early  theory  was  primarily  developed  on  a  qualitative  basis  by 

13  14 
Guinier   '    and  later  extended  to  yield  more  quantitative  results  by 

15 
Guiner  and  Fournet   .   Small  angle  scattering  is  a  phenomenon  commonly 

observed  with  visible  light..   The  halo  observed  around  the  moon  due  to 

scattering  of  the  light  by  small  mist  droplets  is  the  result  of  small 

angle  scattering.  To  observe  such  an  effect,  the  size  of  the  scattering 

particles  must  be  of  the  order  of  ten  to  a  hundred  times  the  size  of  the 

wavelength  of  the  scattered  light.   For  this  type  of  scattering  the 

following  approximate  relationship  is  observed: 

in  which  £  is  the  angle  of  diffraction,  A  the  wavelength,  and  d  the 
mean  diameter  of  the  particles.  The  similarity  between  this  relation 
and  Bragg' s  Law, 

\  -  2d  S»"  6  (3-2) 

is  apparent.   However,  one  must  not  associate  this  type  of  diffraction 


11 


with  that  produced  by  the  regular  periodic  cyrstallites  to  which  Bragg's 
Law  is  applied*   It  is,  in  fact,  the  very  nature  of  the  variation,  or 
non-periodicity  of  structure  which  gives  rise  to  the  small  angle  scatter- 
ing to  which  the  theory  will  be  applied.   Guinier  extended  the  relation- 
ship in  equation  (2- 1 J  to  include  the  same  type  of  diffraction  phenomenon 
in  the  region  of  X-rays.   If  the  wavelength  and  particle  diameter  were 
both  reduced  in  the  same  proportion,  particles  of  the  order  of  ten  to  a 
hundred  angstroms  should  likewise  produce  detectable  small  angle  scatter- 
ing. 

X-rays  are  scattered  by  the  electrons  of  the  diffracting  medium. 
Ordinary  X-ray  patterns  are  obtained  from  the  scattering  of  electrons 
of  atoms  situated  at  regular  coherent  lattice  points.   However,  low- 
angle  scattering  depends  only  upon  the  external  dimensions  of  the  parti- 
cle, that  is,  the  contour  of  the  change  in  election  density  and  not  its 
internal  structure.   The  X-ray  diffraction  of  this  type  ordinarily  ob- 
served is  of  the  form  of  continuous  scattering  at  angle  of  two  degrees 
and  below. 

B.   Theory  -  In  order  to  make  application  of  the  phenomena  to  the  deter- 
mining of  particle  size,  a  relationship  between  scattered  intensity  and 
angle  containing  some  dimensional  parameter  characteristic  of  the  parti- 
cle must  be  developed.   This  parameter  obviously  must  be  sensitive  to 
the  difference  in  the  electron  density  due  to  the  inhomodenity  created 
by  the  particle.   The  derivation  of  appropriate  formulas  will  be  found 
in  Appendix  I,  and  only  the  results  presented  below.   The  characteristic 
particle  dimension  is  designated  by  the  symbol  R.   It  must  be  dependent 
upon  the  electron  distribution  which  in  turn  is  a  function  of  the  atom 
distribution  within  the  particle.   The  total  number  of  scattering  points 
is  therefore  N  atoms  times  the  atomic  number.   The  relative  displacement 

12 


of  these  points  from  some  reference,  in  particular  the  center  of  gravity, 

2 
is  of  interest.    Take  R  as  the  mean  of  the  squares  of  the  distances  of 

each  atom  from  the  center  and  attach  to  it  a  value  equal  to  its  atomic 

weight.   The  value  of  R  then  is  exactly  that  quantity  in  mechanics  call- 

2 
ed  the  radius  of  gyration.   It  is  a  distance  such  that  MR  =  T  where  M  is 

the  particle  mass  and  I  its  moment  of  inertia  relative  to  the  center  of 
gravity.   This  parameter  therefore  is  only  an  indication  of  mass  distri- 
bution and  not  shape.   It  is  because  of  this  very  point  that  the  basic 
theory  still  leads  to  only  qualitative  or  semi-quantitative  results. 
The  majority  of  Guinier's  later  work  "  was  devoted  to  the  refinement  of 
this  aspect. 

Appendix  I  yields  the  relationships: 

for  orientated  particles  ...... 


(3-3) 


Vie  =  N»Z   e^('±A^R^") 

for  randomly  orientated  particles  ...... 

in  which  I  =  scattered  intensity  over  an  angle  £  (in  radians) 

I   =  the  scattering  intensity  by  one  electron 


p 

(obtained  by  Thomson's  formula) 

N  =  the  number  of  particles 

n  =  the  difference  in  the  number  of  electrons  in  the 
particle  and  an  equal  volume  of  the  external 

medium. 

Taking  the  logarithom  of  each  side  of(2-3)one  obtains: 

o 

9         A  TV*  9   9 

Log1QI  a  Log  Nn  Ig 2"   Log1QeR  £ 

A 
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plol  Log      .1  £         shoul  all 


values   of  £     ,    its    slop.-,    m,    being     -  4  TT  /  ,\~    Loq  ..S     |v   .       rhis 
yields   a   value    of 


Rr         AJ2 =r       ,0  76rXfy^ 

for  orientated  particles  and  .418  /\  nT"™   f°r  randomly  orientated 

particles.   Usually  one  plots  Log.  'I  as  a  function  of  the  square  of 

10 

the  magnitude  of  the  vector  h,  where: 

\         4Tr  Sm  6   ^     aire 

•^  I  \  A        for  small  angles 

(9  =    1/2  the  scattering  angle) 

C«   Interpretation  of  results  -   A  straight  line  curve  predicted  by 
equations  i  3-3)  and  (3-4)  is  seldomly  obtained  in  practice.   To  obtain 
a  strictly  linear  curve  demands  only  one  discrete  particle  size  and  of 
course  good  precise  experimental  data.   Figure  3-1  is  a  good  example  of 
these  conditions  for  a  solution  of  albumin.   The  radius  of  gyration  is 
obtained  from  the  slope  of  the  linear  plot.   From  other  information  about 
the  particle  such  as  the  density,  modecular  weight,  and  chemical  struc- 
ture, R  can  be  transformed  into  a  particle  shape  with  specific  dimensions, 
Now  take  the  case  presented  in  Figures  3-2  and  3-3.   What  statements  can 
be  made  about  particle  sizes  for  these  experimental  data?   The  theory  of 
small  angle  scattering  exhibits  its  greatest  validity  for  particles  of 
dilute  concentrations.   This  is  due  to  the  basic  hypothesis  that  it  is 
the  difference  in  electron  density  of  the  particle  and  the  medium  upon 
which  such  scattering  depends,,   For  increasing  concentrations  and  orienta- 
tions a  point  is  reached  in  which  there  appears  to  be  a  degree  of  long 
range  order  and  particle  interference.   This  tends  to  create  a  step  or 
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Fig,  3-1  Low  angle  scattering  of  a  solution  of  albumin  after 
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Fig.  3-2  Particle  interference 


Fig.  3-3  Mixtures  of  particle 
sizes 
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even  a  maxima  in  the  ofcherwi     ntinuous  scattering  curve.   Figure 
3-2  represents  this  situation.   The  most  representative  of  scattering 
patterns  is  represented  in  Figure  3-3.   Only  qualitative  information 
can  be  obtained  from  these  curves.   In  curve  1,  the  variation  of  the  slope 
at  low  angles  is  sufficiently  small  to  permit  an  extrapolation  to  the 
origin.   The  slope  at  the  origin  gives  an  average  radius  of  gyration.   One 
can  say  that  the  sample  does  not  contain  particles  of  too  large  an  R  to 
escape  measurement.   Curves  2  and  3  do  not  permit  such  an  extrapolation. 
The  conclusions  to  be  made  are:   Curve  2,  the  sample  contains  many  large 
particles  and  few  small  ones;  Curve  3,  the.  sample  contains  some  very  large 
particles  but  not  many  small  ones. 

Guinier  '  and  Debye   extend  the  theory  in  an  attempt  to  interpret 
quantitatively  curves  exhibiting  mixtures  of  particle  sizes.   His  analy- 
sis is  based  upon  distribution  functions  in  which  he  applies  Maxwellian 

distributions.   A  more  straight- forward  graphical  procedure  has  been 

17 

developed  by  Jellinek,  Solomon  and  Frankuchen   in  which  the  contribu- 
tions of  particle  size  groups  are  successively  subtracted  out  of  the 
curve.   The  curves  of  Figure  3-3  represent  such  a  sum  of  the  contribu- 
tions of  many  particle  sizes.   They  point  out  that  the  scattering  of  X- 
rays  by  a  homogeneous  assemblage  of  spherical  particles  of  radius  R  can 
be  given  by  the  equation, 

i(e)  =cwi<e  J  5) 

where  I  is  the  scattered  intensity  as  a  function  of  9  (half  the  scat- 
tering angle),  h  =  r— —  ,  W  is  the  weight  of  tie  specimen  doing  the 

A 

scattering  and  C  is  a  constant  for  a  given  material.   The  distribution 
of  particle  sizes  may  be  represented  by  the  function  W(R)  where  W(R)  is 
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the  weight  fraction  within  the  sample  of  radius  R.   Summing  the  total 
scattered  intensity  leads  to  the  expression  for  the  curve: 

x(q)-c)    w(r}r  e  dR  3"6) 

2  2 

The  curve  of  Log  I  (9)  vs,  h  will,  at  the  higher  values  of  h  ,  yield 

a  slope  representing  a  single  value  function  of  R  and  with  an  intercept 

3 
equal  to  CWR  .   Equation  (2-6)  suggests  possibilities  of  obtaining  W(R) 

as  a  continuous  function  of  1(9)  by  application  of  the  Fourier  transform 

18 
method.   A  treatment  of  this  analysis  is  given  by  Bauer   involving  a 

gieat  detailed  process.   The  Jellinek,  Solomon,  and  Frankuchen  method 
greatly  simplifies  matters  by  obtaining  in  place  of  the  continuous  func- 
tion of  W(R)  a  graphically  obtained  set  of  discrete  fractions  of  W(R)  over 
the  observed  range.   Their  data  and  graphical  analysis  for  alumina  gel 
are  presented  in  Figure  3-4s    The  tangent  to  the  curve  at  the  greatest 
observed  angle  of  scattering  is  taken  and  extrapolated  back  to  the 
ordinate  as  the  intercept  K. .   A  point  by  point  subtraction  of  the 
tangent  from  the  curve  yieldsthe  new  curve  (dashed  in  the  figure).   This 
represents  the  scattering  contributed  by  all  particles  of  size  larger 
than  the  value  yielding  slope  m  ,   From  this  curve  the  tangent  m_  is 
taken  and  the  intercept  K?  is  obtained.   The  process  is  repeated  until 
extrapolation  and  point  by  point  curve  subtraction  is  no  longer  feasible. 

The  slopes  are  calculated  and  from  an  appropriate  graph  of  slope  vs.  R 

19 
(Jellinek,  Solomon  and  Frankuchen)    the  discrete  radii  are  obtained. 

3  3 

Since  the  intercepts  are  propcrtional  to  R  W(R),  the  quantity  K./R. 

f"Vi 

represents  the  fraction  of  the  i   particle  contained  within  the  sample. 
The  complete  calculation  for  carbon  black  is  presented  in  Table  1. 
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2   (***) 


Fig.  3-4  Plot  of  Log.  1(G)  against  r   for  Alumina  Gel, 
Jellinek,  Solomon,  and  Frankuchen,  Ind„  Eng. 
Chem.  Annual  Ed.,  18,  172  (1946). 
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Ordinate 

Inte 

rcept 

Ki 

0.61 

K2 

4.3 

K3 

39 

K4 

115 

K5 

140 

K6 

48 

TABLE  I 


Calculation  of  Particle  Sizes  from  Small-Angle  Scattering  Data  for 
Alumina  gel  (Jellinek,  Solomon,  and  Frankuchen?) 


Calculation  of  Slope    0  0         A  r,  ,  ,  .. 

log  I   r* i^Ti   r2  R       R2   ></^  «  ,0  f  ^  • 

'K      Proportion 

0.61   0.484  630  0.000769  10  1,000  6.10  0.194 

1.633  457  0.00358  22  10,600  4.06  0.129 

2.591  240  0.0108  37  50,600  7.71  0.245 

3.061  140  0.0219  54  157,000  7.32  0.232 

3.146  107  0.0294  62  238,000  5.88  0.186 

2.681    33   0.0812    103  1,092,000  .  0.44   0.014 

Total   31.51   1.000 

Jellinek  et  al  in  their  original  publication  acknowledged  the 
apparent  empirical  nature  and  only  piecewise  continuous  size  distri- 
bution of  this  method  of  analysis.   They  could  only  defend  its  validity 

by  stating  that  it  works  for  the  systems  (alumina  gel  and  carbon  black) 

20 
to  which  it  was  applied.   In  a  later  analysis  of  alumina  gel   they 

were  unable  to  obtain  consistent  results. 

21 

Bragg,  Corwin  and  Buttrey   investigated  mixtures  of  silica  sols  of 

two  known  particle  diameters  in  an  attempt  to  verify  methods  of  particle 
size  distribution.   They  were  able  to  substantiate  the  prediction  that 
the  scattering  at  small  angles  arises  predominantly  from  the  larger 
particles  and  that  the  smaller  ones  control  the  scattering  at  larger 
angles.   The  application  of  Jellinek' s  graphical  analysis  was  unsuccessful 
in  reproducing  the  known  particle  sizes  introduced.   They  state,  however, 
that  much  of  their  scattering  data  was  obtained  at  values  of  h  beyond 
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which  the  Guinier  approximation  is  valid.  A  method  due  to  Porod 
which  relates  scattering  at  larger  angles  and  total  surface  area 
was  also  applied.   Porod  predicts  the  following  relationship: 

I(W)  =  ZTC?   Vh4  (3-7) 

in  which  ^  is  the  mean  electron  density,  S  the  total  particle  surface 
area  and  h  as  previously  defined.   In  collimating  systems  of  large 

slit  widths,  the  quantity  h   is  replaced  by  h  .   A  plot  of  In  I  vs  In  h 

3 
should  yield  a  straight  line  of  slope  equal  to  -3  (using  h  )  and  inter- 

2 
cept  of  2n  ^  S.   For  spherical  particle  disperse  system  the  geometri- 
cal relationship  R  -  3  ^  /S        >  where  V  is  Che  particle  volume,  obvi- 
ously holds.   An  analogous  surface  averaged  particle  radius  can  be 
obtained  for  polydispersed  systems.   The  calculated  radius  favors  those 
particles  which  make  the  largest  contribution  to  the  surface  area. 

The  above  survey  substantiates  Guinier' s  original  prediction  that 
the  application  of  small  angle  scattering  in  obtaining  accurate  quantita- 
tive results  is  indeed  limited.   For  widely  separated  spherical  particles 
of  one  diameter,  exact  sizes  can  be  measured.   Non-spherical  particles 
can  be  quantitatively  investigated  if  complementary  information  such  as 
volume,  surface  area,  and  molecular  weight  are  known.   Mixtures  of  sizes 
and/or  unknown  shapes  lead  to  qualitative  results  of  widely  varying 
validity  depending  upon  the  means  of  analysis  employed. 
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CHAPTER  IV 

APPLICATION  OF  SMALL-ANGLE  INFRARED  SCATTERING  TO  PARTICLE  SIZE 

IN  SILICON  CRYSTALS 

A.   General  -   Guinier  developed  his  theory  of  small-angle  scattering 
by  an  extension  of  the  classical  theory  of  atmospheric  scattering  of 
visible  light  down  to  the  region  of  X-rays  and  particles  of  comparable 
size.   It  therefore  seems  reasonable  that  a  similar  extrapolation  can 
be  made  toward  the  region  of  the  infrared  spectrum  and  particle  dimen- 
sions of  the  order  of  microns.   The  limitation  to  materials  which  are 
relatively  transparent  to  infrared  radiation  is  obvious.   An  inspection 
of  equation  (3-1),  £  ~  -~-  suggests  two  possibilities.   One  might  at- 

d 

tempt  to  increase  the  wavelength,  maintaining  a  given  particle  size, 
and  thus  obtain  scattering  at  larger  angles.   This  would  eliminate  the 
complexities  associated  with  small  angle  cameras  and  diffraction  tech- 
niques.  However,  in  the  case  of  X-rays,  longer  wavelengths  are  increas- 
ingly attenuated  by  matter.   This  introduces  problems  associated  with 
reduced  intensities.   The  theoretical  limitation  placed  upon  thir  first 
alternative  is  that  the  expressions  developed  in  the  small-angle  theory 

are  derived  using  an  approximation  called  the  "law  of  Guinier"*   Equa- 

2  2 
-KR  6 
tion  (3-3)  and  (3-4)  for  the  scattered  intensity  contain  the  term  e 

in  which  K  represents  ~r~~  .   As  shown  in  Appendix  I  this  exponential  is 

an  approximation  of  the  series 

_k^c-z  zz     k  **+£.+      k!_r!i8 
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The  approximation  coincides  with  the  exact  expression  only  up  to  the 
4th  power  of  h.   Therefore,  for  large  h,  the  relationship 

is  no  longer  valid.   In  fact,  the  Guinier  analysis  breaks  down  for  hR^.1. 
The  above  approximation  is  the  law  of  Guinier.   Therefore,  one  cannot 
increase  the  wavelength  in  an  attempt  to  obtain  scattering  at  larger 
angles  and  still  apply  Guiner's  small-angle  scattering  theory. 

The  other  suggested  possibility  as  previously  mentioned  is  that  of 
increasing  both  \    and  d,  maintaining  the  same  proportion.   In  this 
manner,  the  angle  limitation  placed  upon  £  is  not  exceeded  and  larger 
size  particles  can  be  measured.   The  study  of  silicon  heterogeneities 

will  be  directed  toward  using  this  relationship. 

22 
Kaiser,  Keck  and  Lange   were  able  to  show  by  means  of  a  9,1  micron 

wavelength  absorption  band  the  existance  of  SiO~  clustering  in  silicon 

crystals.   In  Chapter  II  it  was  pointed  out  that  Schwuttke  was  able  to 

make  some  semi-quantitative  conclusions  of  oxygen  impurity  clustering 

and  dislocation  networks  in  silicon  by  infrared  scattering.   Debye  and 

23 
Lederhandler  "  stated  that  they  were  able  to  obtain  information  on  SiO. 

particle  size  in  silicon  by  measuring  the  angular  distribution  of  infrared 

scattering.   None  of  the  above  investigators  were  able  to  establish  a 

discrete  particle  size. 

The  problem  of  interpretation  of  scattering  data  still  remains  open. 

With  what  kind  of  particle  shape  can  one  associate  dislocations?   Can 

SiCL  clustering  be  considered  to  consist  of  spheroids,  elipsoids  or 

platelets?  There  is  strong  indication  that  oxygen  precipitation  occurs 
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along  dislocation  networks  outlining  a  stringy  network,  or  possibly, 
along  rods  of  an  average  diameter  but  "infinite"  length. 

B.   Dislocations  -  Dash's  decorating  techniques  allows  one  to  draw  some 
conclusions  about  dislocation  networks.   Grown-in  dislocations  tend  to 
form  in  favorable  orientations.   The  formation  of  a  hexagonal  structure 
in  the  (111)  plane  by  the  intersection  of  edge  dislocations  and  screw 
components  was  pointed  out  in  Chapter  II.   Plastic  deformation  induces 
dislocation  loops  from  Frank-Reed  sources  in  the  (111)  plane.   Networks 
of  loops  ranging  from  2  to  50  microns  in  diameter  were  detected.   These 
observations  suggest  the  use  of  equation  (3-3)  for  orientated  particles 
for  low  angle  scattering  analysis.   The  fulfillment  of  the  widely  separ- 
ated particle  criterion  can  be  assured  by  an  inspection  of  the  scattering 
profile.   A  dense  system  will  exhibit  particle  interference  and  a  drop 

off  in  intensity  at  zero  angle.   A  mixture  of  sizes  will  be  indicated  by 

2 
the  concave  upward  nature  of  the  plot  of  Log  I  vs  € 

The  extension  of  small— angle  scattering  theory  to  the  measuring  of 
dislocations  networks  offers  one  gross  possible  shortcoming.  The  scat- 
tering is  due  to  the  difference  in  the  electron  density  of  the  particle 
and  the  surrounding  medium.  Does  a  dislocation  meet  the  description  of 
a  particle  in  this  aspect?  Dislocations  by  definition  do  represent  dis- 
continuities in  atomic  structure.  This  represents  only  a  variation  in 
electronic  density  covering  a  magnitude  of  only  atomic  distances.  How- 
ever, there  is  a  stress  field  surrounding  the  core  which  falls  off  ac- 

2 
cording  to  1/r  ;  r  being  the  distance  from  the  core.   Compared  to  the 

magnitudes  of  atomic  forces,  this  may  be  considered  as  a  relatively 
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long  range  force.   From  the  consideration  of  Hooke's  Law,  this  stress 

24 
is  only  a  manifestation  of  atomic  displacements.   Cottrell    indicates 

that  one-half  of  the  total  strain  energy  surrounding  a  dislocation  lies 
v/ithi;:  a  region  near  the  core.   This  offers  possibilities  of  contribut- 
ing sufficient  variation  in  the  electron  density  to  affect  scattering 
of  the  infrared  radiation.   Furthermore,  Guinier  emphasized  that  "the 
low-angle  scattering  depends  on  the  external  dimension  of  the  particle 
and  not  on  its  internal  structure."   One  can  thus  draw  the  conclusion 
that  a  network  of  dislocations  warrants  the  description  of  a  group  of 
particles. 

As  reported  in  Chapter  II,  Schvuttke  has  reported  the  detection 
of  infrared  scattering  in  silicon  crystals  containing  dislocation 
densities  of  at  least  10  per  cm  .   This  observation  favorably  supports 
the  "particle  nature"  of  dislocation  networks  suggested  above. 

C.   Oxygen  Impurities  -  As  previously  noted  in  Chapter  II,  oxygen  ex- 
ists as  clusters  in  silicon  in  the  form  of  Si0«.   Schwuttke  reported 

17  3 

scattering  in  crystals  containing  at  least  10   oxygen  atoms  per  cm  . 

Kaiser  reported  Rayleigh  scattering  at  concentrations  of  5X10   scatter - 
ing  centers  (10   oxygen  atoms)  per  cm  ,   Both  investigators  made  the 
assumption  that  the  Si09  clusters  are  spherical  in  shape.   This  assump- 
tion has  not  been  verified.   A  great  deal  of  evidence  leads  one  to  be- 
lieve that  precipitation  of  the  oxide  occurs  at  dislocation  sites  as 
do  the  copper  atoms  observed  in  Dash's  decorating  technique.   Micro- 
scopic examination  fails  to  reveal  the  presence  of  a  separate  Si0~ 
phase.   This  however  only  serves  to  insure  that  clusters  greater  than 
approximately  0.1-1  micron  do  exist.   Submicroscopic  spheriodal  particles 
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and  dislocation  site  precipitates  could  still  be  formed.   The  exis- 
tance  of  submicroscopic  spheriodal  clusters  would  obviously  be  most 
favorable  to  the  small  angle  scattering  method.   They  would  tend  to  be 
nearly  uniform  in  size  and  widely  separated  in  distribution.   All  three 
of  these  conditions  (spherical,  uniform,  dilute)  characterize  the 
optimum  conditions  for  the  application  of  the  Guinier  analysis. 

The  alternative  of  the  diffusion  of  oxygen  to  dislocation  sites 

still  exists.   It  is  known  that  dislocation  sites  play  a  significant  role 

25 
in  the  nucleation  of  precipitates.    The  problem  reduces  to:   once 

nucleation  has  commenced,  what  shape  do  the  particles  take?   Continuous 
precipitation  along  dislocation  lines  establishes  the  same  question  of 
meeting  the  particle  criterion  as  discussed  in  section  B. 

The  foregoing  discussion  presents  strong  evidence  of  the  applic- 
ability of  small-angle  infrared  scattering  techniques  to  silicon  crystals, 
With  the  appropriate  complimentary  information,  a  good  estimate  of  parti- 
cle size  and  distribution  should  be  obtained.   Chapter  VI  will  be  devoted 
to  the  analysis  of  experimental  data  by  use  of  the  ideas  presented  here. 
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CHAPTER  V 

EXPERIMENTAL  PROCEDURE 

A.   Specimen  preparation  -  A  1,5  cm.  diameter  silicon  single   rystal 
was  obtained  from  the  Shockley  Transistor  Laboratory.   This  crystal 
was  grown  in  the  ^100/  direction  and  measured  10  ohm- cm  resistivity. 
Several  slices  ranging  from  1.5  to  6  mm.  were  cut  perpendicular  to  the 
growth  axis  and  polished  by  standard  metallographic  procedures  to  minimize 
surface  scratches  and  pitting.   A  check  of  the  finish  at  400  X  magni- 
fication revealed  a  relatively  small  density  of  surface  pitting. 

It  was  initially  anticipated  that  absolutely  parallel  plane  sur- 
faces would  be  necessary  in  order  to  obtain  consistant  results.   Re- 
peated attempts  were  unsuccessful  in  cutting  parallel  faces.   It  was 
decided  that  a  small  amount  of  wedge  could  be  tolerated  if  its  gradient 
is  parallel  to  the  plane  of  detector  rotation.   If  the  taper  were  not  in 
the  plane  of  rotation,  the  refracted  beam  would  contain  a  vertical  com- 
ponent which  would  vary  across  the  face  of  the  radiation  detector.   A 
taper  in  the  plane  of  rotation  merely  caused  a  refraction  which  estab- 
lished the  non-diffracted  or  "zero  angle"  beam,  I  ,  at  a  point  other 
than  directly  in  line  with  the  incident  rays.   Accordingly,  the  crystals 
were  mounted  and  orientated  by  use  of  an  optical  gonimeter.   The  mea- 
sured deviation  from  parallel  averaged  about  1  degree.   This  amount  of 
wedge  was  checked  by  means  of  the  infrared  beam  using  the  relationship 
for  minimum  deviation  of  refracted  light.   For  minimum  deviation: 

S.  i*\  ]fx   d> 
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n  ~  index  of  refraction  of  the  silicon  for  the 
wavelength  used  (3.449  for   A  =  2.3   )  26 

oC =  angle  between  crystal  faces. 

Omin  =  deviation  angle. 

The  refracted  angle,  <£min,  as  measured  by  peaking  intensity,  agreed 

to  within  10'  arc  of  the  value  determined  by  the  optical  goniometer. 

B.   Optics  -  A  photograph  of  the  experimental  set-up  appears  in  figure 
5-1.   The  basic  instrument  is  an  American  Optical  Company  table  spectro- 
meter which  had  previously  been  modified  for  use  as  an  optical  gonimeter. 
In  an  instrument  used  to  measure  small  angle  scattering  there  is  a  pro- 
blem of  achieving  adequate  intensity  without  excessively  sacrificing  re- 
solution.  One  must  make  a  decision  as  to  the  poorest  acceptable  resolu- 
tion and  then  strive  to  obtain  sufficient  intensity  consistant  with  the 
sensitivity  of  the  detector  and  radiation  source.   Two  pairs  of  slits  were 
used;  one  pair  to  collimate  the  incident  beam  and  one  pair  to  shield  the 
detector  from  all  radiation  other  than  the  transmitted  diffracted  beam. 
Slit  separation  was  limited  to  13  cm.  by  the  collimator  barrels  of  the 
spectrometer.   Initially,  the  source  of  radiation  used  was  a  carbon  arc. 
The  carbon  arc  is  a  rich  source  of  near  infrared,  but  difficulty  was  en- 
countered in  maintaining  a  stable,  well  defined  vapor  crater  which 
created  alignment  problems.   A  zirconium  arc  was  then  substituted  for 
the  carbon  arc  with  a  moderate  decrease  in  intensity  but  considerable 
gain  in  stability.   Several  runs  were  made  at  varying  slit  openings  and 
an  acceptable  signal  was  obtained  with  a  width  of  0.79  mm. 

The  wavelength  employed  was  governed  by  the  optical  transmission  of 
silicon  and  the  type  of  detectors  available.   Figures  5-2  and  5-3  represent 
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Fig,    2-2     Transmission  properties  of   silicon,    (^=   15C)IL-cm, 

thickness   0.5  mm  ? 
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the  factors  governing  th.  Lee  ted.  A  wave      h  of  2.3  mi- 

crons was  decided  upon  because  of  the  high  peaking  sensitivity  of  PbS 
photo-conductors  at  this  wavelength.   This  choice  is  consistent  with 
the  broad  transmition  band  of  silicon  in  the  near  infrared  region,   The 
beam  was  monochromated  with  a  "spike"  interference  filter.   As  previous- 
ly mentioned,  obtaining  sufficient  intensity  in  the  beam  was  a  problem. 
Figure  5-4  shows  how  the  total  diffracted  energy  is  successively  dimin- 
ished from  that  of  the  initial  source.   The  energy  radiated  is  represent- 
ed by  the  area  under  the  curve  of  I  =  F  (  /»).   A  detector  responds  to 
total  energy  received  within  the  range  of  its  sensitivity.   The  filter 
greatly  diminishes  the  energy  received  by  the  detector.   Additional  re- 
ductions in  tnergy  are  caused  by  components  in  the  optical  train.   The 
The  filter  transmits  only  approximately  80%  of  the  incident  radiation 
at  the  center  of  the  wave  length.   Another  fraction  of  this  transmitted 
beam  is  absorbed  by  the  crystal. 

An  attempt  was  made  to  increase  the  light  intensity  by  concentrating 
the  beam  with  a  lens  system.   It  was  suspected  that  the  added  attenuation 
due  to  the  glass  optical  elements  would  be  more  than  offset  by  the  beam 
concentration.  Also,  if  the  system  were  successful,  NaCl  or  KBr  lenses 
could  be  substituted  for  glass,  thus  reducing  the  attenuation.  The 
focusing  of  the  light  with  lenses  did  increase  the  detector  response. 
However,  this  gain  was  obtained  at  the  expense  of  collimation  and  associ- 
ated directional  resolution.   One  cannot  increase  the  energy  radiated  per 
area  and  still  maintain  parallel  rays.   In  brief,  the  parallel  rays  can- 
not be  squeezed  closer  together  to  increase  the  photon  density,   To  do  so 
would  appear  to  be  contrary  to  the  second  law  of  thermodynamics.   A  lens 
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system  would  serve  only  in  moving  the  position  of  the  source  as  a  matter 
of  convenience.   The  source  has  a  finite  size  and  cannot  be  concentrat- 
ed.  Therefore,  an  increased  signal  was  obtained  through  increased  slit 
width  and  moving  the  source  as  close  as  practicable,  consistant  with 
adequate  resolution. 

C.   Electronics  -  Associated  with  the  problem  of  the  low  intensity  was 
the  problem  of  signal  amplification.   Figure  5- 5  is  a  schematic  sketch 
of  a  typical  detector  circuit.   In  photoconducting  cells,  it  is  necessary 
to  pass  a  direct  current  through  them  in  order  to  record  the  change  in 
resistance  due  to  photo-conductivity,   Maximum  sensitivity  is  obtained  in 
PbS  cells  with  an  applied  voltage  of  the  order  of  100  volts.   "Dark" 
resistance  of  a  photo-conductor  is  usually  1  megohm.   For  a  maximum 
transfer  of  power,  the  output  or  load  resistor  R  should  be  matched  to 
the  detector  resistance  and  thus  also  1  megohm.   This  creates  a  voltage 
divider  providing  approximately  50  volts  across  the  load  resistor.   With 
the  small  signals  obtained,  of  the  order  of  millivolts,  it  was  not  possible 
to  separate  the  variation  in  voltage  from  the  voltage  measured.   Increas- 
ing the  load  resistor  or  using  d.c.  amplifiers  merely  increased  the  V  a- 
long  with   V  maintaining  the  same  ratic   Direct  current  amplifiers  have 
been  used  with  PbS  cells  with  large  signals  but  not  of  the  order  of  magni- 
tude encountered  in  the  experimental  setup.   Separation  of  the  detector 
signal  from  the  steady  d.c.  bias  can  be  achieved  by  use  of  chopped  radia- 
tion as  indicated  in  Figure  5  6   An  a.c  amplifier  is  used  to  amplify 
the  a.c.  component  created  by  the  chopper.   Amplifier  units  were  obtained 
from  a  Perkins-Elmer  infrared  spectrometer.   Essentially,  all  the  compon- 
ents except  for  the  optical  train  were  easily  adapted  for  use  in  an  in- 
frared diffraction  system.   The  units  consist  of  a  13  cycle  per  second 
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chopper,  pre  !  second  amplifier  with  rectified 

output,  and  a  Leeds  and  Northrup  recorder.   With  proper  shielding  of 
the  detector,  a  high  signal  to  noise  ratio  was  obtained  and  relatively 
low  intensity  signals  adequately  amplified. 
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CHAPTER  VI 
EXPERIMENTAL  RESULTS 

Scattering  data  were  obtained  for  three  slices  of  the  original 

crystal.   Thicknesses  were  2,  3,  and  6  mm.   Dislocation  densities 

7   8       2 
obtained  by  etch-pit  counts  were  estimated  to  be  10  -10  per  cm   for 

each  specimen.   Since  the  first  tests  were  conducted  without  heat 
treating  the  crystals,  no  contribution  to  the  light  scattering  by  oxy- 
gen was  anticipated,,   Microscopic  examination  did  not  reveal  the 
presence  of  a  separate  phase  due  to  Si09  clustering  or  other  impurities. 
Figures  6-1  and  6-2  represent  the  plotted  scattering  data  for  the 

2  mm  and  6  mm  specimens.   The  3  mm  crystal  produced  essentially  identical 

2 
data  as  the  2  mm  specimen.   The  highest  value  of  €   correspond  to  a 

scattering  angle  of  2  degrees.   Also  indicated  on  each  curve  is  the  limit 
specifying  Rh  iS  1.   Therefore  the  data  was  obtained  at  angles  well  with- 
in those  for  which  the  Guinier  analysis  is  applicable.   A  check  was  made 
for  parasitic  diffraction  due  to  the  collimating  slits.   This  was  per- 
formed with  the  optical  arrangement  identical  to  that  of  the  test  runs 

but  without  the  crystal.   The  intensity  was  reduced  to  essentially  zero 

2 
at  angles  of  1  degree  (  €  —4).   Sharp  minima  were  observed  at  an  angle 

of  40  minutes.   This  caused  the  discontinuity  in  the  specimen  scattering 

2 
curves  observed  at  values  of  6.   equal  to  2-3.   The  main  conclusion  drawn 

from  the  parasitic  diffraction  test  is  the  absence  of  any  contribution  to 

2 
scattering  in  the  region  of  interest,  £  >4,  due  to  the  slits. 

From  the  limiting  slopes,  radii  of  gyration  of  5.70^  and  4.79yU 
are  calculated  for  the  2mm  and  6  mm  specimens  respectively.   To  inter- 
pret these  values  in  terms  of  a  particle  size  auxiliary  information  of 
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particle  shape  is  required.   As  pointed  out  in  Chapter  II,  the  results 
of  earlier  investigators, suggest  a  pol)  ^onal  or  circular  platelet 
configuration.   The  difference  in  electron  density  is  concentrated  about 
the  periphery  of  the  loop.   A  good  approximation  of  particle  shape  would 
be  that  of  an  annul us  with  a  radius  of  gyration  of: 


R 


If  the  width  of  the  dislocation  line  and  bordering  displaced  atoms  is 
considered  to  be  narrow  compared  to  the  radius  of  the  loops  one  can 
ssume  r  O^.    r_.   For  this  approximation,  the  expression  reduces  to: 


as? 


r  =  R 
o 

In  order  to  calculate  a  particle  size,  we  must  therefore  assume  that  the 
particles  are  discs  of  negligible  thickness  whose  radius  of  gyration  is 
the  effective  radius  in  the  plane  of  the  disc.   This  assumption  leads  to 
a  particle  size  with  diameters  of  11.4^and  9.58^<for  the  specimens  of 
figures  6-1  and  6-2  respectively.   These  figures  are  of  the  same  order 
of  magnitude  of  the  values  observed  by  Dash  using  decoration  techniques. 

The  resistivity  measurement  imply  a  considerable  impurity  content.   The 

2 
presence  of  oxygen  will  be  seen  at  large  values  of  £  .   The  results  in- 
dicate no  particles  less  than  11  L\    in  diameter.   We  therefore  conclude 
that  the  impurities  are  in  solution  in  the  silicon  crystal  or  at  most 
trapped  in  the  dislocations. 
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Figures  6-1  and  tte  a  distribution  of  particle  sizes. 

The  long  linear  portion  of  the  slope  at  large  angles  describe  a  larger  ,;er 
centage  of  smaller  particles.   The  sharp  bend  in  the  curve  at  small  angles 
shown  in  figure  6-1  indicate  a  smaller  proportion  of  large  particles  for 
the  2  mm  specimen.   However,  the  data  at  angles  less  than  1  degree  are 
not  completely  reliable  due  to  the  contribution  of  parasitic  scattering 
by  the  collimating  slits* 
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APPENDIX  I 

DERIVATION  OF  SMALL-ANGLE  SCATTERING  FORMULAS 

A.   Scattering  By  Several  Waves 

A  sinusoidal  wave,  of  frequency  oO,  is  defined  by  its  amplitude 
A  and  phase  0   : 

a  —   A  Cos  C  <^>~t  +  &) 

This  can  be  expressed  in  complex  variable  notation: 

in  which  a  is  the  real  part  of  U.  .   For  several  waves  of  same  frequency 
propagated  in  the  same  direction  the  observed  amplitude  will  be  the  sum 
of  the  amplitudes  of  each  individual  wave. 

U=Eu;=    ejut2A.ejA' 

The  resultant  amplitude  and  phase  can  be  obtained  by  vector  addition 
from  the  following  diagram: 


F,<v  *-l 


Since  the  intensity,  I,  is  equal  to  the  square  of  the  amplitude,  it  is 
evident  from  the  above  diagram  that : 

I  -    (?  A,   Cos  0A  f  +    (?A<  Sf*  (f>A) 
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Expressed  in  complex  variable  notation 

I  =  0  /.  u*    =     (f 


e^2,A,  e'-fl  e'J^2kAk  e." 


=r   5.  2k  A,  Ak  e.^-*")  =  ?;  2k  A,  AK  Cos (<t>A  -tf)     i -i 


B.   Phase  Difference  of  Waves  Scattered  by  Two  Electrons 


A  g 


F.v  I"2 


In  the  above  figure,  S   represents  a  unit  vector  along  the  incident  rays. 
Points  0  and  A  represent  two  electrons  scattering  the  incident  rays.   S 
is  the  scattered  ray  unit  vector.   The  phase  difference  between  the  rays 
scattered  at  0  and  A  can  be  represented  by  the  retardation  distance  MA  +•  An, 
From  vector  algebra: 


|mAl     =     OA  •  S, 


1A 


f\ 


— 5  — 4> 


l«y»A  *  An[        =       OA  •   (  S0  -  S   ) 


The  phase  difference  is  therefore: 


0    = 


2TT 

A 


OA 


.  ( s„-s) 


Since  |s   -  Sjis  equal  to  £Sin  Q   and  is  parallel  to  ON  « 


4V 


0  -  ^-r    a  b'  i  V\  9  Co  s  m 

A 


(X-2) 
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C .   Intensity  of  Low-Angle  Scattering  by  Small  Particles 

1.   Oriented  particles.   The  intensity  of  the  scattering  by  an 
individual  atom  is  I  f  where  I   is  the  intensity  of  the  scattering  by  one 
electron  and  f  the  atomic  scattering  factor.   The  amplitude  is  then  re- 
presented by   A  ~  v  Prt  £   and  formula  (I-l)can  be  written; 

■o 


P13.  T-3       C 
In  figure  1-3,   £   represents  the  scattering  angle  and  is  equal  to  20. 
Since  the  scattering  considered  here  is  at  very  small  angles,  OD  may  be 
regarded  as  perpendicular  to  S  for  all  values  of   c.   .   The  length  a 
is  the  projection  of  A  on  ON.   Therefore  the  phase  difference,,  formula 

K. 

(1-2),  can  be  written: 

0k  =  -  £E  £  .  Sak  Cos  ja    =:     ^  ^  ^  Oak 

and  .        „^ \ 

r  -  xt  1 1  fj,  C*  cos  (  *f  €  •  aK  a,) . 

Letting  0.^0.^  A^  ^   and   -a-L.  £  —  y\       ,  the  intensity  can  be  expressed 


as  a  cosine  series: 
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*/^«&tf- t***tef«<fck)  - 


The  various  scattering  points,  a  ,  a a  can  be  projected  on  the 


normal  OD  as  X,,  X-...X  .   Hence 
1    2    p 
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"* 'A  K  •  ^ K 


The  double  summation  5   s?    $     ( ,  (    jP  \  *         !  first  over  k  is 

\  %  fJAuV"  =  £  2*  *;  x,<L  -  *  £  K>  2k  *fc  **   +  ft  tf  *  fh 

Choose  the  origin  such  that   £  f   x .   ^  Q  ,      This  represents 

the  projection  of  the  center  of  gravity  of  the  particle  on  0D„   The 
second  summation  gives: 

<£   \v    is  the  total  number  of  electrons,  n,  in  the  particle,  and 


the  length  R.  is  defined  such  that: 


Since  f   is  approximately  proportional  to  the  mass  of  the  atom, 

2 
where  M  is  the  total  mass  of  the  particle  and  MR   the  moment  of  inertia 

of  the  particle  with  respect  to  the  plane  normal  to  OD  and  passing 

through  the  center  of  gravity.   The  double  summation  is  then  written: 


and 


where 


r/r    r  n^[  i  -   h1^  +  - 


] 


2X. 


X 

This  is  the  expansion  of  the  exponential  series  e  .   Thus  the  forrrula 
for  the  scattering  of  an  isolated  particle  is:        z.       •)        -> 

which  is  formula  (3-3)  in  Chapter  III, 

2.   Random  oriented  particles  -  For  randomly  oriented  and  irregular- 
ly placed  particles,  the  total  intensity  is  the  sum  of  the  intensities 
diffracted  by  the  individual  particles.   When  there  are  N  particles: 


7rer  ^n[i-^ 


N 


+ 
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2 
where  MR  '  etc.  are  the  moments  of  inertia   with  respect  to  the  three 


where    r- —       represents  the  mean  of  the  squares  of  the  values  of 

R  .   The  mean  value  of  the  moment  of  inertia  of  a  particle  in  relation 

to  planes  of  random  orientation  passing  through  its  center  of  gravity 
must  be  found.   From  classical  mechanics: 

Ml<    ^'(m^   ^Vm^)   +>1/m*£) 

.'i 

principle  planes  and   c^ ,  ^  .  V*   the  direction  cosines  to  the  normal 

2 
of  the  plane  considered.   The  mean  value  of  R.   is: 

2 
MR  being  the  moment  of  inertia  in  relation  to  the  center  of  gravity, 

R  is  the  radius  of  gyration.   Therefore: 

which  is  the  series  expansion  of   .  •>      ■}  "L    -*   •, 

r/rer  n5N  e.     3     ^   nlN  e     *a^k  € 

This  is  formula(3-4)  of  Chapter  III. 
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